function P_cross = func_LogitCross(Response,X,X0,n,p,lambda_BCFVH,K)

    P_hat = zeros(n,1); 

    for k=1:K    
        Psi_phi = max(eps,sqrt(mean((repmat((Response-P_hat),[1,p]).*X).^2,1)));
        [B,FitInfo] = lassoglm(X/diag(Psi_phi),Response,'binomial','Link','logit','Lambda',lambda_BCFVH/n,'Standardize',0);                    
        P_hat = 1-1./(1+exp(FitInfo.Intercept+(X/diag(Psi_phi))*B));            
    end
                    
    P_cross = 1-1./(1+exp(FitInfo.Intercept+(X0/diag(Psi_phi))*B)); 

end